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Carbon nanoscroll is a unique topologically open configuration of graphene nanoribbon possessing 
outstanding properties and application perspectives due to its morphology. However molecular 
dynamics study of nanoscrolls with more than a few coils is limited by computational power. Here, 
we propose a simple model of the molecular chain moving in the plane, allowing to describe the 
folded and rolled packaging of long graphene nanoribbons. The model is used to describe a set of 
possible stationary states and the low-frequency oscillation modes of isolated single-layer nanoribbon 
scrolls as the function of the nanoribbon length. Possible conformational changes of scrolls due to 
thermal fluctuations are analyzed and their thermal stability is examined. Using the full-atomic 
model, frequency spectrum of thermal vibrations is calculated for the scroll and compared to that 
of the flat nanoribbon. It is shown that the density of phonon states of the scroll differs from the 
one of the flat nanoribbon only in the low (cu < 100 cm“^) and high {to > 1450 cm“^) frequency 
ranges. Finally, the linear thermal expansion coefficient for the scroll outer radius is calculated from 
the long-term dynamics with the help of the developed planar chain model. The scrolls demonstrate 
anomalously high coefficient of thermal expansion and this property can find new applications. 


I. INTRODUCTION 

During last decades, various carbon nanostructures 
have attracted an increasing attention of researchers, 
due to their unique electronic, mechanical, and chem¬ 
ical properties, as well as many potential applications. 
Numerous studies of single-layer graphene sheets and 
graphene nanoribbons (GNR) have started to take place 
in recent years [1-5]. Secondary graphene structures like 
folds or scrolls can be placed in a separate class of carbon 
nanomaterials whose existence is ensured by the action of 
relatively weak van der Waals bonds between sp^-bonded 
monatomic carbon layers. The observation of scrolled 
graphite plates under surface rubbing was first reported 
in 1960 [6]. The authors of this work have suggested that 
the lubricating properties of graphite are due essentially 
to the rolling up of packets of layers, which then act like 
roller bearings provided a low coefficient of friction. The 
thickness of the scrolls was estimated as of order of 100 
planes. 

The spiral shape and geometric parameters of carbon 
nanoscrolls (CNS) are determined by the balance of en¬ 
ergy gain due to increase of the number of atoms involved 
in van der Waals interactions with the energy loss due to 
graphene bending. Several experimental techniques for 
obtaining and study of CNS have been reported. [7- 
14] Properties of CNS have also been studied in a se¬ 
ries of theoretical investigations. Electrical, optical and 
mechanical properties of short CNS have been described 
from ah initio calculations. [15-17] Mechanical properties 
of CNS and various scenarios of their self-assembly have 
been described by means of molecular dynamics method. 
[18-29] 

Mechanical properties and the lowest vibration fre¬ 
quency of long CNS have been described in the frame¬ 
work of the continuum model of a spiral elastic rod[19, 
24, 30, 31], where the bending energy of the rod is com¬ 


pensated by the energy gain from the interaction of ad¬ 
joining walls. 

The resonant oscillation of a CNS near its fundamental 
frequency might be useful for molecular loading/release 
in gene and drug delivery systems. [19] Internal hollow 
core is one of the main structural features of carbon nano¬ 
scrolls. Due to this cavity the system of coils at low tem¬ 
peratures can serve as an effective storage of hydrogen 
atoms, [32, 33] whereas a separate scroll can be used as 
an ion channel. [34] 

Under lateral compression hollow cored configuration 
of the scroll demonstrates a weak resistance and the core 
can collapse at a moderate load. This feature opens 
the perspective for application of parallel-stacked CNS 
as an efficient device sensitive to pressure, which may 
be used as nano-sized pumps and filters. [31, 35] Simula¬ 
tions buckling and post-critical behavior of CNS under 
axial compression, torsion, and bending has revealed the 
occurrence of kinks and folds. [36] 

The molecular dynamics technique has proved to be 
a very powerful tool for simulation of mechanical prop¬ 
erties and deformations mechanisms of CNS. However, 
the full-atomic models are very demanding in computa¬ 
tional power making consideration of long-term dynamics 
of CNS with a large number of coils almost impossible. 

Addressing these challenges, in this paper we propose 
a simple model of the planar molecular chain capable of 
describing the longitudinal and flexural motion of CNR 
and allowing the study of folded and/or scrolled CNR. 

In Sec. II the chain model of the carbon nanoribbon is 
introduced and the parameters of the model are fitted to 
some results in frame of the full-atomic model. In Sec. 
Ill the chain model is applied to simulate the secondary 
structures of single layer nanoribbon, such as folded 
and scrolled configurations. Then graphene nanoribbon 
scrolls are analyzed in more details in Sec. IV. Erequency 
spectra of the flat nanoribbon and nanoribbon scroll are 
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Figure 1: A scheme for constructing the chain model of 
the carbon nanoribbon, (a) Full-atomic model of GNR with 
zigzag orientation; (b) the chain model. Atomic rows of the 
nanoribbon oriented along the 2 -axis are numbered by the 
index n. 


calculated in Sec. V using the full-atomic model. Ther¬ 
mal expansion of scrolls is analyzed in Sec. VI using the 
chain model. Section VII concludes the paper. 


II. CHAIN MODEL OF THE CARBON 
NANORIBBON 

Graphene nanoribbon is a narrow, straight-edged 
stripe of graphene. It is well-known that graphene is elas¬ 
tically isotropic material and thus, its longitudinal and 
flexural rigidity depend weakly on chirality. For definite¬ 
ness, GNR with the zigzag orientation will be considered 
as shown in Fig. 1 (a). 

In the flat configuration the nanoribbon is supposed 
to lay in the iz-plane of the three-dimensional space. 
The nanoribbon can be described as a periodic structure 
with the step a = roCOs(7r/6), where tq = 0.1418 nm is 
the C-C equilibrium valence bond length. Let us con¬ 
sider such modes of the nanoribbon motion in which the 
carbon atoms in the atomic rows parallel to the 2 -axis 
move as the rigid units only in the xy-plane. Under this 
assumption, tensile and flexural GNR dynamics can be 
described by the chain of point-wise particles moving in 
the a;?/-plane. Atomic rows of the nanoribbon oriented 
along the 2 -axis are numbered by the index n as shown 
in Fig. 1 (b). 

The chain model of GNR schematically shown in Fig. 2 
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Figure 2: Chain of particles numbered by the index n in the 
ly-plane, modeling tensile and bending motion of graphene 
nanoribbon. Potentials V and U describe the longitudinal and 
bending stiffness of the chain, respectively. The potential W 
(not shown here) describes the weak van der Waals bonds 
acting between layers of the chain in the folded or rolled con¬ 
formations. 

can be described by the following Hamiltonian 

N N-1 

n—1 n—1 

N-1 N-4 N 

+ + E (1) 

n—2 n—1 k—n-\-4 

where ?/„ are the coordinates of n-th particle, r„ = 
|v„| is the distance between particles n and n -I- 1, with 
vector v„ = {xn+i — Xn, yn+i — Un) Connecting these par¬ 
ticles, On is the angle between vectors v„ and —Vji_i, Vnk 
is the distance between particles n and k (index i = 1 if 
the difference — n is an odd number and i = 2 ii k — n 
is an even number). 

The Hrst term in Eq. (1) gives the kinetic energy of 
the chain with M = I2mp being the mass of carbon atom 
{nip = 1.6603 ■ 10“^^ kg is the proton mass) and the 
dot denotes differentiation with respect to time t. The 
harmonic potential 

V{r) = ^K{r-af, (2) 

describes the longitudinal stiffness of the chain. 

The angular anharmonic potential in Eq. (1) 

C/(0) = e[cos(0)-I-1], (3) 

stands for the flexural rigidity of the chain with cos(0„) = 
— (v„_i, v„)/r„_i7’„ being the cosine of the n-th ’’valent” 
angle. 

The potential Wi{rnk), (* = 1,2) in Eq. (1) de¬ 
scribes the weak van der Waals interactions between 
particles n and k, located at the distance Xnk = 
\/{xk — XnY + {Vk — VnY- These interactions, acting 
between nanoribbon layers, must be taken into account to 
describe folded or scrolled conformations of GNR. With¬ 
out these interactions the only stable configuration of a 
nanoribbon is the flat one. 

Parameters of the Hamiltonian K and e are determined 
to fit the dispersion curves of the flat carbon nanoribbon 
presented in Eig. 1(a) to the dispersion curves of the 
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straight chain model shown in Fig. 2. To do so, we 
consider dynamics of the flat nanoribbon, in line with 
the chain model, assuming that the atoms can move only 
in the xy-plane with all atoms in the rows along 2 -axis 
displaced equally. Let us denote the coordinates of the 
n-th atomic row as Xn and ?/„ and introduce the notation 
for the vector u„ = (xmljn)- For the flat nanoribbon the 
weak van der Waals interactions do not contribute to the 
dynamics and can be neglected. Then the Hamiltonian 
for the nanoribbon can be written in the form 


+ CXD 


E 


n— — oo 


M(u„,u„) -h P(u„_i,u„,u„+i,u„+ 2 ) 


( 4 ) 

The first term in Eq. (4) describes the kinetic energy, 
while the second one stands for the potential energy of 
interatomic interactions, both per one carbon atom of n- 
th atomic row, located far from the nanoribbon edges. 
Thus, the effect of nanoribbon edges is not taken into 
account or, in other words, the nanoribbon width effect 
is not taken into account. 

To describe the carbon-carbon valence interactions let 
us use a standard set of molecular dynamics potentials 
[37]. The valence bond between two neighboring carbon 
atoms a and (3 can be described by the Morse potential 


Ui{ua,up) = ei{exp[-ao(r- ro)] - 1 }^, r = ju^ - u^j], 

( 5 ) 

where ei = 4.9632 eV is the valence bond energy and 
ro = 0.1418 nm is the equilibrium valence bond length. 
Valence angle deformation energy between three adja¬ 
cent carbon atoms a, /3, and 7 can be described by the 
potential 

U2{Ua,Up,U^) = e2(c0SV3 - COS(/Jo)^, (6) 


where cos(p = (u^ - u^, - U/jj/du^, - u^j • |u^ - u^l), 

and (/ 5 o = 27r/3 is the equilibrium valent angle. Param¬ 
eters oq = 17.889 nm“^ and £2 = 1.3143 eV can be 
found from the small amplitude oscillations spectrum of 
the graphene sheet [38] . Valence bonds between four ad¬ 
jacent carbon atoms a, j3, 7 , and 5 constitute torsion 
angles, the potential energy of which can be defined as 


Uzicj)) = €^{1-cos (j)), (7) 


where (j) is the corresponding torsion angle ((/> = 0 is the 
equilibrium value of the angle) and £3 = 0.499 eV. 

A detailed discussion of the choice of the interatomic 
potential parameters can be found in [37]. The same set 
of potentials has been successfully used to simulate the 
heat transfer along the carbon nanotubes and nanorib¬ 
bons [39, 40] for the analysis of spatially localized oscilla¬ 
tions [41-46] and also for the investigation of theoretical 
strength and post-critical behavior of deformed graphene 
[47-50]. 

Hamiltonian Eq. (4) generates the following set of the 
equations of motion 


“t“P2(Un—1; ll^_|_2) -f F3(ll^_2, 11^, Uyj^i) 

-fP 4 (u„_ 3 ,U„_ 2 ,U„_i,U„), ( 8 ) 

where vector-function 
d 

Fk = ^ —P(U1,U2,U3,U4), fc = l,2,3,4. 
duk 

It is convenient to use the relative coordinates of atoms 
w„(t) = u„(t) — u°, where u° are the equilibrium coor¬ 
dinates. For the analysis of small-amplitude vibrations 
(|w„| <g; ro) we use the following linearized equations of 
motion 


- Mw„ = BiWn + P2(w„_i -H W„+i) 

+B3{Wn-2 + W„+ 2 ) -f P4(W„_3 -f- Wf„+ 3 ), (9) 


where matrices Pi = Pn -I- P 22 + P 33 + P 44 , P 2 = P 12 + 
P 23 + ^ 34 , P 3 = Fi 3 + P 24 , P 4 = Pi4, and matrix 


Fki 


d 

dv/kdwi 


P(0, 0,0,0). 


We seek for the solution of linear system Eq. (9) in the 
form of the wave 


w„(t) = Aexp[i(gn — wt)], (10) 

where oj is the frequency of the wave, A is the ampli¬ 
tude vector, q € [ 0 , 7 r] is the dimensionless wavenumber. 
Substituting Eq. (10) into the linear system Eq. (9) we 
obtain the dispersion relation 

jPi -I- 2 cos(q)B 2 + 2 cos{2q)B3 -|- 2 cos( 3 g)P 4 — w^P| = 0, 

( 11 ) 

where E is the unity matrix. 

The dispersion relation Eq. (11) is the second order 
polynomial with respect to the squared frequency 
The corresponding dispersion relation has two branches 
0 < ujy{q) < uJx{q) as plotted in Fig. 3 by the solid lines. 
The low-frequency branch w = ujy(q) describes the dis¬ 
persion of the transverse plane waves when lattice nodes 
leave the nanoribbon plane and move along y axis (bend¬ 
ing nanoribbon vibrations). The high-frequency branch 
w = uJz{q) represents the dispersion of the longitudinal 
plane waves when the nodes move along a:-axis (longitu¬ 
dinal nanoribbon vibrations). 

Velocity of long-wavelength plane phonons coresponds 
to the velocity of sound which is 

Vx = a liin u!x(q)/q = 17510 m/s 

g ->0 

for the longitudinal phonons and 

Vy = a lim ujy{q)/q = 0 
q —>0 

for bending phonons. 

Similarly, we can obtain the dispersion curves for the 
chain model. In this case, potential energy in the Hamil¬ 
tonian Eq. (4) is defined as 


Miln = Fi(Un,Un+l,Un+2,Un+3) 


P (Un—1 ; Un 1 ^n+l) — ^(^n) “1“ ^(^n)- 
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Figure 3: Dispersion curves of the wide carbon nanoribbon. 
Curve 1 stands for the longitudinal (tension-compression) and 
curve 2 for the transversal (bending) phonons. Solid lines 
correspond to the full-atomic model, while dashed lines to 
the chain model with appropriately chosen parameters. 

Potential parameters ii'of Eq. (2)andeofEq. (3) should 
be chosen in a way to achieve the best fit of the disper¬ 
sion curves obtained for the full-atomic model. We are 
interested in the long-wavelength modes of the nanorib¬ 
bon motion and thus, the best fit should be achieved in 
the range of small values of q. The choice 

K = 405 N/m, e = 3.50 eV, (12) 

assures the coincidence of the longitudinal and flexural 
rigidity of the chain model and the nanoribbon, as can 
be seen in Fig. 3, where the result for the chain model is 
plotted by the dashed lines. 

The van der Waals interactions between carbon atoms 
are described by the pairwise Lennard-Jones potential 

Wo{r) = 4e[(CT/r)^^ - {a/rf], (13) 

where r is the distance between two carbon atoms. 
The parameters of the Lennard-Jones potential e = 
0.002757 eV and a = 0.3393 nm were fitted to repro¬ 
duce the interlayer binding energy [51], interlayer spac¬ 
ing [52, 53], and c-axis compressibility [54] of graphite. 
The potential Eq. (13) with these parameters is shown 
in Fig. 4 by the curve 1. 



Figure 4: The pairwise Lennard-Jones potential Wo{r) de¬ 
scribing the van der Waals interaction between two carbon 
atoms (curve 1). Solid lines 2 and 3 show the potentials Wiir) 
and W^ir) (see in the text), while their approximation by the 
modified Lennard-Jones potentials Eq. (16) are shown by the 
dashed lines 2 and 3. 

The long range interaction between the chain nodes 
n and k is described by the van der Waals interactions 
of the atoms belonging to n-th and fc-th atomic row of 
the nanoribbon. Therefore, the interaction energy can be 
expressed as 

OO 

Wi(r) = J2 l^o{rj,i) + Wo{r,, 2 )], (14) 

j^-OC 

+ (-0-5?'o + 

rj,2 = P + (-l-5ro -h 3jro)^]^/^, 
if the difference fc — n is an odd number and 

OO 

W2{r) = ^ [lEo(r,-3) + Wo(r,-4)], (15) 

j=-oc 

rj,3 = [r^ + (3Jro)^]^/^ 

^1.4 = + {-ro + 3jro)^]^/^ 

if fc — n is an even number. 

Interaction potentials Eq. (14) and Eq. (15) are well 
approximated by the modified Lennard-Jones potential 

Wi(r) = 4e,[(cr,//(r))i2 - {cJi/f{r)f], 

f{r) = r,{rln)°‘\ r* = (16) 

For z = 1, the modihed potential parameters are £i = 
0.008652 eV, ui = 0.31636 nm, and ai = 0.86. For i = 2, 
parameters are £2 = 0.008029 eV, 0^2 = 0.32607 nm, 
and 02 = 0.90. The interaction potentials of nanorib¬ 
bon atomic rows Eq. (14) and Eq. (15) together with the 
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corresponding modified Lennard-Jones approximations 
Eq. (16) are shown in Fig. 4 by the solid and dashed lines, 
respectively. Practically perfect coincidence of these po¬ 
tentials can be observed. 

Summing up, for GNR shown in Fig. 1 (a), we have 
developed the chain model depicted in Fig. 2 and de¬ 
scribed by the Hamiltonian Eq. (1) with the parameters 
fitted to reproduce the long-wavelength phonon spectrum 
of GNR (see Fig. 3) and the van der Waals interactions 
acting between carbon atoms in folded or scrolled confor¬ 
mations of the nanoribbon (see Fig. 4). The chain model 
describes only such modes of GNR deformation in which 
the atomic rows parallel to z-axis move as rigid units only 
in the xy-plane but not in ^-direction. The nanoribbon 
width effect is not taken into account. 



III. SECONDARY STRUCTURES OF SINGLE 
LAYER NANORIBBON 

To find stable structures of a one-layer carbon nanorib¬ 
bon the following minimization problem should be con¬ 
sidered 


N-l N-1 

e = y . + E 

n—l n—2 

N-A N 

+ E E Wi{rnk) -t min : (17) 

n=l fc=ri+4 


where the minimization of potential energy of the chain 
model having N nodes is performed with respect to the 
coordinates of the nodes u„ = (xmUn), n = 1,2, 

The nanoribbon length is defined by the number of nodes, 
L = {N -l)a. 

The energy minimization was carried out numerically 
using the conjugate gradient method. In order to check 
the stability of the resulting stationary configuration 
we calculate the eigenvalues of the 2N x 2N 
matrix of the second derivatives 


B = 


dE 


dUn,idUk,j 


N, 


N, 


f,,0 \N j 

J m = l/ n=l,Z = l,fc=l,J = l 


(18) 


The stationary chain configuration is stable only if all 
eigenvalues of a symmetric matrix B are non-negative: 
Ai > 0, i = 1,2,...,2Y. Note that for stable configura¬ 
tion the first three eigenvalues are always zero Ai = A 2 = 
A 3 = 0. These eigenvalues correspond to the rigid mo¬ 
tion of the chain in the ly-plane, with two translational 
and one rotational degrees of freedom. The remaining 
positive eigenvalues Ai > 0 correspo nd to the oscillation 
eigenmodes with frequencies Wi = ^J\i/M, i = 4,..., 2N. 

Stationary structure of the chain depends on the ini¬ 
tial configuration used to solve the minimization prob¬ 
lem Eq. (17). Changing the initial configuration, a 
variety of stable configurations can be found. Linear 


Figure 5: Stable stationary conformations of the chain model 
with N = 140 nodes representing the nanoribbon of length 
L = 17.070 nm; (a) rolled, (b) double-folded, (c) triple-folded, 
and (d) rolled-collapsed. Not shown here is the straight stable 
configuration (flat nanoribbon). 


chain configuration, representing flat nanoribbon, is al¬ 
ways stable. The weak van der Waals interactions be¬ 
tween nodes give rise to the existence of other, more ad¬ 
vantageous in energy, stationary states of the chain in 
the two-dimensional space. As exemplified in Fig. 5, the 
chain consisting of Y = 140 nodes and having length 
L = {N — l)a = 17.070 nm, in addition to the flat state, 
can be stable in (a) rolled, (b) double-folded, (c) triple- 
folded, and (d) rolled-collapsed conformations. Potential 
energy per node, Eq = E/N, is used to compare the en¬ 
ergy of different chain conformations. In case of Y = 140, 
the flat structure has Eq = —0.00453 eV. For the other 
forms one has: Eq = —0.01395 eV for the rolled state, 
—0.01214 eV for the double-folded, —0.00352 eV for the 
triple-folded, and —0.00662 eV for the rolled-collapsed 
state. To understand these figures, one should keep in 
mind that formation of van der Waals bonds lowers the 
structure total energy, while the large curvature regions 
increase the energy. The rolled packing is the most en¬ 
ergetically favorable among the studied conformations of 
the nanoribbon. All the studied non-flat structures have 
energy lower than the flat one, except for the triple-folded 
one. This is explained by the fact that the triple-folded 
structure possesses two loops with large curvature having 
no van der Waals bonds, and such loops have relatively 
large energy. 

The dependence of the normalized energy Eq for dif¬ 
ferent stationary nanoribbon packages on its length L 
is shown in Fig. 6 . In the range L < 5.77 nm the 
planar structure is the only stable configuration of the 
nanoribbon. For L > 5.77 nm, stable rolled structures 
exist. Ghains with L > 6.02 nm {L > 10.19 nm) can 
support stable double-folded (triple-folded) structures. 
The rolled-collapsed structure requires the chain length 
L > 10.19 nm. 
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Figure 6: Potential energy per node, Eq = E/N, as the 
function of the chain length L for flat nanoribbon (curve 
1), double-folded (curve 2), triple-folded (curve 3), rolled- 
collapsed (curve 4), and for the rolled structures with one-coil 
(curve 5), two-coil (curve 6), three coil (curve 7), four coil 
(curve 8) and more coiled (curve 9) structures of the nanorib¬ 
bon. 


The flat nanoribbon has the lowest energy for L < 
10.93 nm. For the nanoribbon length in the range 
10.93 < L < 13.39 nm the lowest energy is observed for 
the double-folded configuration and for L > 13.39 nm 
the most energetically favorable is the rolled structure 
(nanoribbon scroll). 


IV. GRAPHENE NANORIBBON SCROLLS 

In the preceding Section it was shown that GNR hav¬ 
ing length L > 13.39 nm have the lowest energy in the 
rolled conformation among the other studied configura¬ 
tions. That is why, here we focus on the study of nanorib¬ 
bon scrolls. The cross-sectional view of the minimum en¬ 
ergy scroll structures for nanoribbons of increasing length 
can be seen in Fig. 7. The scroll cross-section appears in 
the form of the truncated Archimedes spiral always hav¬ 
ing inner cavity. The scroll structure is determined by the 
balance of energy gain caused by increasing the number 
of atoms having van der Waals bonds with the others 
and the energy loss due to the increase of nanoribbon 
curvature. 

The center of mass can be considered as the center of 



Figure 7: Examples of the stable, equilibrium scroll struc¬ 
tures of the nanoribbon: (a) single coil scroll (nanoribbon 
length L — 7.00 nm, energy per node Eq — 0.00881 eV, num¬ 
ber of coils Nc = 0.96, inner and outer radii Ri = R2 = 
1.11 nm); (b) double coil scroll {L = 7.25 nm, Eq = 0.0123 eV, 
Nc = 1.19, Ri — 0.97 nm, R2 = 1.021 nm); (c) double coil 
scroll {L = 19.53 nm, Eq = —0.01684 eV, Nc = 1.98, Ri = 
1.40 nm, R 2 = 1.73 nm); (d) three coil scroll (L = 19.53 nm, 
Eo = -0.01692 eV, Nc = 2.19, Ri = 1.24 nm, R2 = 1.63 nm); 
(e) three coil scroll {L = 34.26 nm, Eq = —0.02815 eV, 
Nc = 2.96, Ri = 1.51 nm, R2 — 2.17 nm); (f) four coil scroll 
(L = 34.26 nm, Eo = -0.02816 eV, Nc = 3.14, Ri = 1.38 nm, 
R2 = 2.10 nm); (g) four coil scroll (L = 51.45 nm, Eq = 
-0.03432 eV, Nc = 3.94, Ri = 1.58 nm, R 2 = 2.57 nm); (h) 
five coil scroll {L = 51.45 nm, Eo = —0.03429 eV, Nc = 4.05, 
Ri = 1.51 nm, R 2 — 2.54 nm). 

the scroll: 

1 ^ 

"0 = vE"- 

n—l 

where u° = (a;°,j/°) is the two-dimensional radius- 
vectors of the n-th chain node of the equilibrium scroll. 
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In the polar coordinate system it can be written as 

= xo + i?° cos(<;i°), = yo + sin(^!)° ), (19) 

where (x° — a:o)^ + (y° — yo)^ and the discrete 

angle monotonously increases with increasing node 
number n = The spiral can be characterized 

by the number of coils 


- ?5'i)/27r. 


It is also convenient to define the integer number of coils 
m = [Nc] + 1, where [cc] is the integer part of x. Let us 
define the inner radius of scroll by its first coil: 


Ri 


1 

ni 




where ni is the maximal value of index n wherein (j)n < 
+ 27r. The outer radius of the scroll can be defined by 
its last coil as: 


i?2 


1 


N 


iV - 712 + 1 


i—N—n2 


where n 2 is the minimal value of n where (j)n > (j^N — Stt. 

The twisting rigidity of the scroll is characterized by 
the lowest natural frequency wi = y/X 4 /M. This fre¬ 
quency corresponds to the periodic twisting/untwisting 
oscillations of the scroll. In the approximation of a con¬ 
tinuous elastic rod this oscillation motion has been stud¬ 
ied in [19, 30]. 

Let us describe how the scroll structure and the low¬ 
est natural frequency depend on the chain length (see 
Fig. 7 and Fig. 8). Various conformations can be natu¬ 
rally characterized by the number of coils Nc-, energy per 
node Eq, inner and outer radii i?i, i? 2 . 

Single-coil configuration is the only possible in the 
chain length range 5.772 < L < 7.000 nm [Fig. 7 (a)]. 
Double coil configuration is stable in case of 7.245 < 
L < 20.508 nm [Fig. 7 (b), (c)]. The chain length 
range 18.052 < L < 36.718 nm corresponds to sta¬ 
ble three coil scrolls [Fig. 7 (d), (e)]. For the case of 
33.771 < L < 51.542 nm the four-coiled scrolls are ob¬ 
served [Fig. 7 (f), (g)]. If nanoribbon length L > 51.0 nm, 
the scrolls with five or more coils exist [Fig. 7 (h)]. 

One can see that for nanoribbons with the same length 
two stable configurations are possible. For example, in 
the length range 18.052 < L < 20.508 nm stable two or 
three coil scrolls exist [Fig. 7 (c) and (d)]. Nanoribbons 
of length 33.771 < L < 36.718 nm can be packed in three 
or four coil scrolls [Fig. 7 (e) and (f)]. The reason for 
bistability is the result of the interaction of the nanorib¬ 
bon ends. One stable configuration is when the ends are 
close to each other and another one is for somewhat over¬ 
lapped ends. The degree of overlapping decreases with 
increasing nanoribbon length, as it can be seen in Fig. 7. 
Increase of the chain length leads to weakening of the 



Figure 8: (a) Effect of the nanoribbon length L on (a) the 

number of coils of the stationary scroll Nc, (b) inner Ri (lower 
curve) and outer R 2 (upper curve) radii of the scroll, (c) the 
lowest natural vibration frequencies of the scroll, cJi and oJ 2 - 
The curves 1, 2, and 3 correspond to a single-coil scroll; curves 
4-7 to the two-coil scroll; 8-11 to the three-coil scroll; 12-15 
to four-coil scroll; curves 16-19 to the case of five and more 
layered structure of the nanoribbon scroll. The dashed lines 
are the fitting curves: (a) Nc = 0.33Z/°'®®, (b) Ri — 0.67, 
R 2 ~ 0.4L° '^^, and (c) toi = 06/L, with L, Ri, and R 2 given 
in nanometers. 


interaction between ends and thus to weakening of the 
bistability of the scroll packing. 

Increase in the nanoribbon length L results in the 
monotonous increase in the number of coils Nc according 
to the power law Nc Ri 0.33L°®^, see Fig. 8 (a). The 
inner scroll radius Ri increases much slower with L than 








the outer radius i? 2 : Ri ~ 0.67L°-^, R 2 « 0.4L'’-^’^, see 
Fig. 8 (b). Here L, i?i, and R 2 are given in nanometers. 

The eigenmode having lowest positive frequency wi is 
the twisting-untwisting mode when the atoms move along 
the Archimedes spiral. The second and the third low¬ 
est eigenfrequencies correspond to lateral compression- 
extension of the scroll. The scroll symmetry is lowered 
by the nanoribbon ends and for this reason the lateral 
compression in the two orthogonal directions is charac¬ 
terized by the close (but not equal) frequencies W 2 and 
wa. These frequencies depend on L nonmonotonically, 
see Fig. 8 (c). For wi the general trend is the reduc¬ 
tion of the frequency with the growth in L according to 
the law uii ss 36/L for L —>■ 00 . This is in line with 
the asymptotic behavior obtained analytically in Refs. 
[19, 30]. 



Figure 9: (a) Two dimensional chain model and (b) corre¬ 

sponding three dimensional full-atomic model of the nanorib¬ 
bon scroll with the nanoribbon length = 36.84 nm and 
width Lj = 2.55 nm {N = 301, K = Q, number of atoms 
Nm = 2NK = 3612). 


V. FREQUENCY SPECTRUM OF THE 
NANORIBBON AND NANORIBBON SCROLL 


Let us perform the full-atomic three-dimensional mod¬ 
eling of the dynamics of the nanoribbon scrolls to verify 
the two-dimensional chain model. 

Let the set of two dimensional vectors {u° = 

1 ) 2 )}^=! is the solution of the minimization prob¬ 

lem Eq. (17), describing a scroll packing of the nanorib¬ 
bon of length = (iV—l)a. For the nanoribbon of width 
Lz = SKtq (the translational cell consist of 4AT carbon 
atoms), then the coordinates of the atoms in n-th cell are 


^n,4(fc —1) + 1 
■^n,4(/c— 
Xn,i{k-l)+2 
Zn,'l(k-1) + 1 
^n,4(fc—1)+3 
Zn,i(k-1) + 1 
^n,4(fc—1)+4 
2n,4(fc-l)-|-l 

k 


'^ 2 n-l,li yn,i{k-l) + l — W 2 „_i_ 2 i 
3(fc - l)ro, 

'^ 271 , 1 ; yn,iik-l)+2 = '^ 2^,21 

ro/2 -I- 3(fc - l)ro, 

^2n,li 2/n,4(fc-l)-|-3 = ^2^,21 (20) 

3ro/2 + 3(fc - l)ro, 

'*^2n-l,li 2/n,4(fc-l)-|-4 = "^271-1,2) 

2ro -b 3(fc - l)ro, 

1,2,..., A, n = l,2,...,7V/2, 


where, as above, tq is equilibrium C-C valence bond 
length. In the case of odd number of cells, N, the 
nanoribbon consists of Nm = 2NK carbon atoms. The 
two dimensional chain model and the full-atomic nanorib¬ 
bon scroll are presented in Fig. 9 for the nanoribbon 
length = 36.84 nm and width = 2.55 nm {N = 
301, K = 6, number of atoms Naii = 2NK = 3612). 

A set of interaction potentials (5), (6), (7), (13) was 
used for modeling of the nanoribbon dynamics. Valence 
bonds between neighboring atoms in the graphene plane 
are described by the Morse potential (5), valence and 
torsional angles by the potentials (6) and (7). Weak van 
der Waals interactions between scroll coils are described 
by the Leonard-Jones potential (13). Let us consider 
the edge carbon atoms chemically modified by hydrogen 


atoms and thus, having mass of Mq = ISm^, under the 
assumption that the interaction potentials for edge and 
internal atoms are the same. 

Dynamics of the nanoribbon having size = {N — 
l)a X SiFro is described by the Langevin equations 

dH 

Afn,iU„ i = —--rM„ iU„ i-I-i, (21) 

OUn,l 

n = l,2,...,iV, / = l,2,...,4iF, 

where = {xn,i,yn,h Zn^i) is the three dimensional 
radius-vector of (n,/)-th atom, ; is the atom mass 
(Mn^i = I 2 mp for the internal atoms and Mnj = 13mp 
for the edge atoms). Here H is the nanoribbon Hamilto¬ 
nian, F = 1 /tr is the friction coefficient, (velocity re¬ 
laxation time is tr = 0.4 ps), random forces vectors 

are normalized as follows 

(^n,/,i(H)C7n,fcj (I 2 )) — 2,AIkBT5nm^lk^ij^{tl ^ 2 )- 

The set of equations of motion Eq. (21) is integrated 
numerically. Initial conditions corresponding to the sta¬ 
tionary scroll packing of the two-dimensional chain model 
Eq. (20) were used. We took K = 6 and N = lOI, 201, 
301, and 401 to simulate the nanoribbon of width = 
= 2.5524 nm and length = {N — l)a = 12.28, 
24.56, 36.84, 49.12 nm. Numerical simulation at low tem¬ 
perature r = 3 K has shown the absence of any notice¬ 
able changes in the initial structure of the scroll, i.e. the 
stationary configuration of the three dimensional scroll 
is well described by the equilibrium configuration of the 
two-dimensional chain. This confirms the high accuracy 
of the two-dimensional chain model. 

In order to find the phonon density of states for the 
full-atomic flat nanoribbon and for the full-atomic scroll, 
the Langevin equations Eq. (21) were integrated for 10 ps 
to achieve the state of thermal equilibrium at the desired 
temperature, and then the thermostat was switched off 
and free dynamics of atoms was studied. It was demon¬ 
strated numerically that the flat and scrolled nanorib¬ 
bons of length Lx = 36.84 nm are both stable in the 
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Figure 10: Phonon density of states calculated from 
thermal vibrations of atoms for the flat nanoribbon (red 
curves 1 and 3) and for the scroll (blue curves 2 and 4) at 
temperatures (a) T — 300K and (b) T = 900K. The nanorib¬ 
bon length and width are = 36.84 nm and = 2.55 nm, 
respectively. 


temperature range 0 < T < 900 K. The use of full- 
atomic model allows us to find the time dependence of the 
particle velocity on time and then the density of 

phonon states p(a;) normalized such that p{uj)du} = 1. 
The density of phonon states was determined from 600 
homogenously distributed atoms and 256 independent re¬ 
alizations of the initial thermalized nanoribbon state in 
order to increase the calculation accuracy. The result 
for the nanoribbon of length = 36.84 nm and width 
Lz = 2.55 nm is shown in Fig. 10 for the flat nanoribbon 
(red curves 1 and 3) and for the scroll (blue curves 2 and 
4) at temperatures (a) T = 300K and (b) T = 900K. As 
one can see, the frequency spectra of the flat and scrolled 
nanoribbons are very close. Certain difference can be ob¬ 
served only in the low uj < 150 cm“^ and high cj > 1450 
cm“^ frequency intervals. In the range w < 50 cm“^ the 
scroll has phonon density more than two times smaller 
than the flat nanoribbon. This is due to the fact that 
the rigidity of scroll is higher than that of nanoribbon 
and the low-frequency bending and torsional vibration 
modes are absent in the scroll. At high frequencies a 
small blue shift (by 5 cm~^) of the oscillation frequencies 
is observed for the scroll. Scrolling of the nanoribbon 
leads to a moderate increase of oscillations frequencies 
in the range w > 1450 cm“^ due to the van der Waals 
interactions of atoms belonging to adjacent layers of the 
scroll. 


VI. THERMAL EXPANSION OF SCROLLS 

The scroll structure is stabilized by the weak van der 
Waals bonds acting between coils. Thermal fluctuations 
weaken such bonding leading to partial untwisting or full 
opening of the scroll. Fully opened scroll transforms to 
the flat GNR, while partial untwisting results in a reduc¬ 
tion of the number of coils and in a growth of the scroll 
diameter. 

The full-atomic model does not allow to simulate 
the long-term dynamics of wide multi-coiled nanoribbon 
scrolls due to the computer capacity limitations. For 
example, thermal expansion of the scroll can hardly be 
treated by the full-atomic model and this problem is ad¬ 
dressed here in frame of the two-dimensional chain model. 

For the simulation of thermal vibrations of the chain 
the Langevin equations were used 

dH 

Mu„ =--rMu„-bS„, (22) 

dUn 

n = 1,2,..., V, 

where u„ = (xn,yn) is the radius-vector of n-th node, 
H is the Hamiltonian of the chain Eq. (1), N is the 
number of nodes in the chain, F = 1/A is the friction 
coefficient (velocity relaxation time is A = 5 ps), and 
= (Cn,i,Cii, 2 ) is the two dimensional vector of nor¬ 
mally distributed random forces, normalized as 

(fny(A)Cmj(A)) = 2MkBTSnniSijS{ti - A) 

(here ks is the Boltzmann constant). 

The set of equations of motion Eq. (22) was integrated 
numerically. The stationary state of the scroll was used 
as an initial configuration. 

Thermal stability of the scroll depends on the nanorib¬ 
bon length L. The longer is the nanoribbon the larger 
is the energy of van der Waals bonds per atom and the 
higher is the thermostability. From simulations, it was 
found that single-coiled scroll of the nanoribbon having 
L = {N — l)a = 12.28 nm (number of nodes N = 101) 
is stable only for T < 330K and at higher temperatures 
it fully opens in less than t = 10 ns. The two-coil scroll 
with L = 24.56 nm {N = 201) was found to be stable 
within the whole studied temperature range T < 960K. 
Only partial untwisting is observed for the scrolls with 
this and higher values of L in this temperature range. 

Let us take the scroll of the nanoribbon of length L = 
110.40 nm (number of nodes N = 900) for the study of its 
dynamics in the range of temperature 30 < T < 960K. 
The dependence of number of coils A/ and inner and 
outer radii i?i, i ?2 of the scroll on time is shown in Fig. 11. 
As one can see from the graph, thermal vibrations lead 
to decrease of the number of coils Nc and growth of the 
radii i?i and i? 2 - 

Typical scroll configurations at different temperatures 
are shown in Fig. 12. At elevated temperatures the spiral 
structure is maintained and the inner and outer radii of 
the scroll increase with temperature. 
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Figure 11: Time evolution of (a) number of coils Nc, (b) 
outer radius R 2 , and (c) inner radius Ri for the scroll of the 
nanoribbon having length L = 110.40 nm, evaluated at the 
temperatures T = 30 K (curves 1, 2, and 3), T = 300 K 
(curves 4, 5, and 6), and T = 900 K (curves 7, 8, and 9). 



Figure 12: Typical scroll configurations for the nanoribbon 
having length L = 110.40 nm (number of chain nodes N= 
900) at temperatures (a) T — 30K, (b) T = 300K, and (c) 
T = 900K. 


In about 1 ns thermalization of scroll is complete and 
it obtains its equilibrium configuration. Further integra¬ 
tion of the equations of atomic motion allows to deter¬ 
mine the averaged values of the coil number and inner 
and outer radii of the scroll, Ri and i? 2 , correspond¬ 
ing to the given temperature. This values, as the func¬ 
tions of temperature, are plotted in Fig. 13. It can be 
seen that the number of coils decreases linearly, while 
the radii demonstrate a linear increase with tempera¬ 
ture as Nc{T) Ri 6.52 - 0.00I42T, Ri{T) « 7 ?° -b CiT, 



Figure 13: The temperature dependence of (a) the average 
number of coils Nc, (b) the outer radius R 2 , and (c) the inner 
radius Ri of the scroll of the nanoribbon having length L — 
110.40 nm. Dotted lines present the linear approximations. 


where i = 1,2, = 1.782 nm, i?® = 3.63 nm, ci = 

0.00075 nm/K, and C 2 = 0.00047 nm/K. The relative in¬ 
crease in the outer radius of the scroll of the nanoribbon 
having length L = 110.40 nm is R 2 {T)/R 2 w 1 -b cT, 
where the coefficient of linear thermal expansion is equal 
to c = C 2 /i ?2 = 1-3 ■ 10“^K“^. It was found that c 
depends on L such that c is higher for smaller L. For ex¬ 
ample, for L = 73.56 nm one has c = 1.4 • 10“^K“^; for 
L = 36.84 nm c = 1.8 • 10“^K“^; and for L = 25.56 nm 
c = 3.7- lO-^K-i. 

Note that the coefficient of linear thermal expansion 
calculated for the graphene scroll outer radius is two 
(one) orders of magnitude larger than that for graphite in 
a (c) direction (see [55] and references therein reporting 
on the experimental data) and two orders of magnitude 
larger than that for diamond [56]. 


VII. CONCLUSIONS 

In this paper, the two-dimensional chain model (see 
Fig. 2) was developed to accurately and effectively de¬ 
scribe the dynamics of folded and rolled conformations 
of graphene nanoribbons. The Hamiltonian of the model 
is given by Eq. (1) and it takes into account the tensile 
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and bending rigidity of the nanoribbon, as well as the van 
der Waals interactions between layers of the nanoribbon. 
The chain model describes only such modes of nanorib¬ 
bon deformation [see Fig. 1 (a)], where the atomic rows 
parallel to z-axis move as rigid units only in the xy-plane 
but not in z-direction. The nanoribbon width effect is not 
taken into account. Parameters of the chain model were 
fitted to reproduce the low-frequency part of the phonon 
dispersion curves of the flat graphene nanoribbon (see 
Fig. 3). The van der Waals interactions were fitted by 
the modified Lennard-Jones potentials (see Fig. 4). 

The validity of the chain model was demonstrated by 
comparison of the structure of the stationary nanoribbon 
scrolls with the results of full-atomic simulations. 

Potential energy per atom was calculated for flat, 
rolled, double-folded, triple-folded, and rolled-collapsed 
conformations of nanoribbon as the function of its length 
L (see Fig. 5 and Fig. 6). It was found that in the range 
L < 5.77 nm the planar structure is the only stable con¬ 
figuration of the nanoribbon. For L > 5.77 nm, sta¬ 
ble rolled structures exist. Chains with L > 6.02 nm 
(L > 10.19 nm) can support stable double-folded (triple- 
folded) structures. The rolled-collapsed structure re¬ 
quires the chain length L > 10.19 nm. The flat nanorib¬ 
bon has the lowest energy for L < 10.93 nm. For the 
nanoribbon length in the range 10.93 < L < 13.39 nm 
the lowest energy is observed for the double-folded config¬ 
uration, and for L > 13.39 nm the most energetically fa¬ 
vorable is the rolled structure (nanoribbon scroll). These 
results can be easily interpreted in terms of the competi¬ 
tion between energy release due to the formation of van 
der Waals bonds and energy absorption due to bending 
of the nanoribbon. 

Since GNR having length L > 13.39 nm have the low¬ 
est energy in the rolled conformation among the other 
studied configurations, the nanoribbon scrolls were stud¬ 
ied in detail (see Fig. 7 and Fig. 8). It was found 
that single-coil configuration is the only possible in the 
chain length range 5.772 < L < 7.000 nm. Double 
coil configuration is stable for 7.245 < L < 20.508 nm. 
The chain length range 18.052 < L < 36.718 nm cor¬ 
responds to stable three coil scrolls. For the case of 
33.771 < L < 51.542 nm the four-coiled scrolls are ob¬ 
served. For nanoribbon length L > 51.0 nm, the scrolls 
with five or more coils exist. Increase in the nanoribbon 
length L results in the monotonous increase in the num¬ 
ber of coils Nc according to the power law Nc ~ 0.33L° ®^. 
The inner scroll radius Ri increases much slower with L 
than the outer radius R 2 . i?i ^ 0.67i°-^, i ?2 ~ 0.4L° '^^. 
Here L, i?i, and i ?2 are given in nanometers. The 
twisting-untwisting eigenmode having lowest positive fre¬ 
quency LOi was calculated as the function of nanoribbon 
length. For long nanoribbons the asymptotic law was 


found uji « 36/L for L —>■ oo, which is in line with the 
earlier theoretical studies [19, 30]. 

Full-atomic model was used to calculate phonon den¬ 
sity of states for flat and scrolled nanoribbons of length 
Lx = 36.84 nm and width = 2.55 nm at 300 and 900 K 
(see Fig. 10). It was shown that the phonon spectra for 
the two conformations are very close in the entire fre¬ 
quency range. Small difference can be observed only in 
the low oj < 150 cm“^ and high oj > 1450 cm“^ frequency 
intervals. This can be explained by the higher rigidity of 
scroll due to formation of the van der Waals bonds be¬ 
tween adjacent layers of the scroll. The low-frequency 
bending and torsional vibration modes are absent in the 
scroll. 

One of the most important findings of the present 
study has emerged from the application of the developed 
chain model to the simulation of the long-term dynamics 
of nanoribbon scrolls at different temperatures. It was 
found that the relative increase in the outer radius of the 
scroll of the nanoribbon having length L = 110.40 nm is 
characterised by the coefficient of linear thermal expan¬ 
sion of c = 1.3 • 10“'^K“^. It was found that c depends 
on L such that c is higher for smaller L. For exam¬ 
ple, for L = 73.56 nm one has c = 1.4 • 10“^K“^; for 
L = 36.84 nm c = 1.8 • 10“^K“^; and for L = 25.56 nm 
c = 3.7 • 10“"^K“^. The coefficient of linear thermal ex¬ 
pansion calculated for the graphene scroll outer radius 
is two (one) orders of magnitude larger than that for 
graphite in a (c) direction [55] and two orders of magni¬ 
tude larger than that for diamond [56] . Such anomaly in 
the coefficient of thermal expansion can be used in the 
design of nanosensors or other nanodevices. 

The developed planar chain model can help to address 
problems related to the dynamics of open graphene struc¬ 
tures not treatable by the full-atomic simulations. The 
results obtained in frame of the chain model could pro¬ 
vide a better understanding of the mechanical properties 
of CNS-based nanodevices. 
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